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摘 要 : 本 文通 过 SSTk-o 潮流 模型 和 ZGB 空 化 模型 对 某 种 涡轮 泵 进行 了 空 化 流 场 计 算 ， 得 到 了 其 水 
力 性 能 曲线 和 空 化 性 能 曲线 ， 对 设计 工 况 下 涡轮 泵 内 部 流 场 和 压力 脉动 特性 进行 了 分 析 。 结 果 表 明 ， 涡 
轮 泵 易 在 三 个 位 置 发 生 空 化 ,诱导 轮 前 后 的 压力 脉动 主要 由 诱导 轮 叶片 通过 频率 决定 ， 蜗 壳 内 部 压力 脉 
动 受 动静 干涉 影响 较 大 。 
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Fluctuation at Rated Condition in High-speed Turbopump 
Zhang Hao-Chen Zuo Zhi-Gang Liu Shu-Hong Wu Yu-Lin 

(State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China) 
Abstract: This paper used SST k-o turbulence model and ZGB cavitation model to simulate 
cavitation flow field in a turbopump. Its hydraulic performance curves and cavitation 
performance curves were obtained and the inner flow distribution and pressure fluctuation 
characteristics in this turbopump were analyzed under the design point. Results show that there 
are three cavitation regions, and the pressure fluctuation in the front of inducer and behind the 
inducer is determined by the inducer blade passing frequency. Rotor-stator interaction impacts 


the pressure fluctuation in the volute. | . . . | 
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0 引言 

对 于 液体 火箭 发 动机 ， 涡 轮 泵 是 关键 部 件 之 一 ， 但 是 随 着 涡轮 泵 向 高 速 大 功率 密度 发 展 ， 
对 于 涡轮 泵 的 空 化 性 能 有 着 更 高 的 要 求 。 在 涡轮 泵 中 ， 诱 导 轮 是 提高 主 泵 空 化 性 能 的 主 
要 途径 ， 虽 然 诱导 轮 可 以 在 一 定 空 化 条 件 下 运行 ， 但 是 仍然 会 产生 稳定 性 问题 。1999 年 
本 H-I 火箭 发 射 器 发 射 失败 ， 发 动机 残骸 在 海中 发 现 ， 根 据 分 析 ， LE-7 发 动机 的 涡轮 
泵 诱导 轮 中 发 生 了 空 化 不 稳定 性 现象 ， 导 致 发 动机 出 现 故 障 ， 最 后 火箭 发 射 器 坠落 患 ; 
P&W 公司 为 NASA 研发 航天 飞机 主 发 动机 的 改进 型 高 压 液 氧 涡轮 泵 , 研制 中 空 化 引起 了 
严重 的 振动 现象 ， 导 致 诱导 轮 叶片 被 破坏 ,试验 也 因此 提前 停止 站。 因此 研究 涡轮 泵 内 部 
空 化 流 场 及 压力 脉动 有 助 于 其 内 部 流动 稳定 性 的 分 析 和 优化 。 
由 于 涡轮 泵 空 化 实验 不 易 开 展 ， 因 此 目前 研究 其 内 部 流动 的 方法 主要 是 数值 模拟 。M. ML 
Athavale 等 人 口 使 用 了 CFD-ACE+ 软 件 ， 采 用 全 空 化 模型 和 标准 满 流 模型 k-s 模型 ， 其 中 
蒸汽 输 运 方程 是 基于 Rayleigh-Plesset 方程 和 其 它 条 件 导出 ， 计 算 了 SSME 涡轮 和 泵 的 诱导 
轮 和 离心 泵 的 流 场 。 该 文通 过 捕捉 流 场 中 细微 变化 ， 发 现 了 在 不 同情 况 下 空 泡 长 度 的 变 
化 和 空 泡 所 在 的 区 域 。Hosangadi A 等 人 由 和 Dorney D 等 人 采用 标准 kg 模型 对 单一 诱 
导 轮 进行 了 数值 计算 ， 模拟 了 其 内 部 空 化 ,分 析 其 不 稳定 问题 原因 。 本 文通 过 使 用 Fluent 
软件 对 某 一 涡轮 泵 内 空 化 流 场 进行 了 模拟 ,分析 了 其 内 部 空 化 区 域 及 产生 原因 ， 并 且 讨 
论 了 影响 内 部 压力 脉动 的 原因 。 
1 数值 计算 方法 

本 文 对 滴 轮 泵 进行 建 模 ， 将 其 划分 为 进口 段 、 诱 导 轮 、 离 心 轮 、 蜗 过、 出 口 段 五 部 分 ， 
全 流 道 几何 模型 见 图 1。 采 用 六 面体 对 进出 口 及 蜗 壳 进行 网 格 划分 ;采用 四 面体 对 诱导 轮 
和 离心 轮 进行 网 格 划 分 ， 对 于 诱导 轮 叶 尖 和 管 壁 之 间 的 间隙 进行 加 密 。 涡 轮 泵 主要 参数 
及 计算 边界 条 件 见 表 1。 为 了 保证 结果 的 准确 性 ， 本 文 进行 了 网 格 无 关 性 验证 ， 见 图 2， 
本 文 使 用 547 万 网 格 进行 计算 。 

表 1 涡轮 泵 参数 及 计算 条 件 


Tab.1 Parameters of turbopump and simulation conditions 
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涡轮 泵 参数 
设计 转速 (r/min) 40000 
设计 流量 (L/s) 3.1 
进口 压力 (Pa) 202650 
诱导 轮 叶 片 数 2 
离心 轮 叶 片 数 6 长 6 短 
计算 条 件 
相 流 模型 混合 物 模 型 
汕 流 模型 SST k-o 模型 
空 化 模型 ZGB 模型 
进口 条 件 速度 进口 
出 口 条 件 压力 出 口 
压力 速度 耦合 SIMPLEC 算法 
差分 方法 二 阶 迎 风格 式 
温度 20°C 
工 质 水 /水 蒸汽 
时 间 步 长 (s) 2.083x10* 
收敛 残 差 10” 
进口 诱导 轮 


蜗 壳 离心 轮 下 
\ 挫 500 
400 600 
出 口 网 格 数 /万 个 


图 1 涡轮 泵 全 流 道 几何 模型 ” 图 2 网 格 无 关 性 验证 

Fig.1 Profile ofturbopump ~ Fig.2 Grid-independent verified 

2 定常 计算 结果 分 析 

本 文采 用 定常 数值 模拟 方法 ， 预 测 了 涡轮 泵 的 外 特性 曲线 ， 如 图 3 所 示 。 对 于 本 文 的 涡 
轮 泵 来 说 ， 在 约 1.1Q0 附近 处 ， 效 率 出 现 最 大 值 ; 在 设计 工 况 下 ， 效 率 为 67%， 设 计 工 
况 下 的 扬程 达到 了 532m。 本 文 改变 进口 压力 ， 得 到 涡轮 泵 的 空 化 性 能 曲线 (图 4)。 临 界 
空 化 点 一 般 取 效率 下 降 3% 的 工 况 点 ， 从 图 中 可 以 发 现 设计 工 况 对 应 的 是 临界 空 化 点 ， 推 
测 本 文 使 用 的 涡轮 硝 仍 存在 优化 的 空间 。 
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诱导 轮 两 个 叶片 的 吸力 面 和 压力 面 的 压力 分 布 如 图 5 所 示 ， 可 以 看 出 诱导 轮 内 的 低压 
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有 两 个 位 置 ， 一 个 是 诱导 轮 叶 片 进口 轮 缘 处 ， 男 一 个 是 诱导 轮 叶 片 出 口 靠近 轮 载 处 。 可 


以 从 图 中 看 出 ， 诱 导 轮 两 个 叶片 的 吸力 面 的 压力 变化 趋势 是 一 致 的 ， 压 力 面 亦 如 此 ， 同 
时 压力 面 的 压力 比 吸力 面 大。 


c 
图 5 诱导 轮 压力 分 布 


Fig.5 Pressure in the inducer blades 
图 6 显示 诱导 轮轴 截面 的 流 线 分 布 情况 ， 可 以 看 出 在 诱导 轮 叶 片 进口 处 有 回流 涡 存 在 ， 
这 是 诱导 轮 的 进口 回流 现象 。 诱 导 轮 的 进口 回流 包含 了 垂直 轴 面 方向 的 流 涡 和 绕 流 线 的 
回旋 流 这 两 类 。 这 是 因为 诱导 轮 旋 转 的 叶片 对 流体 的 作用 力 不 均 匀 ， 由 于 离心 力 的 作用 ， 
在 诱导 轮 叶 片 外 缘 处 与 轮 慌 处 之 间 会 产生 压 差 ， 但 是 由 于 叶片 始终 旋转 着 ， 这 种 运动 压 
差 也 一 直 保 持 着 ， 因 此 会 造成 流体 从 叶片 进口 外 缘 倒 流 到 进口 管 段 ， 这 些 回 流 部 分 在 主 
流 中 会 使 得 主流 做 旋转 运动 。 同 时 在 实际 运行 状态 下 ， 实 际 液 流 角 会 比 设计 工 况 下 的 液 
流 角 小 一 些 ， 这 会 使 得 诱导 轮 叶片 进口 前 缘 产 生 不 均匀 的 圆周 速度 分 量 ， 圆 周 速度 分 量 
产生 了 绕 流 线 的 流 涡 。 诱 导 轮 进口 回流 会 造成 局 部 低压 。 
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图 6 诱导 轮回 流 涡 

Fig.6 Backflow vortex in inducer 

设计 工 况 下 涡轮 泵 流 道 内 部 的 空 化 区 域 见 图 7a。 空 化 区 域 主要 发 生 在 三 个 地 方 ， 分 别 在 
诱导 轮 叶片 入 口 轮 缘 处 ， 诱 导 轮 出 口 靠 近 轮 载 处 ， 离 心 轮 长 叶 厂 吸力 面 进口 处 ， 这 三 个 
区 域 都 是 低压 区 。 在 相同 工 况 下 ， 计 算 了 无 诱导 轮 时 离心 叶轮 内 部 的 空 化 流 场 。 对 于 离 
心 叶轮 内 空 化 区 域 〈 图 7p) ， 可 以 看 到 没有 诱导 轮 的 条 件 下 ， 离 心 叶轮 内 部 空 化 从 离心 
轮 进口 一 直 延 伸 到 各 个 流 道 的 中 后 段 。 通 过 对 比 可 以 发 现 ， 诱 导 轮 能 改善 离心 叶轮 内 的 


空 化 性 能 。 


图 7 涡轮 泵 空 化 区 域 
Fig.7 Cavitation regions in turbopump 
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分 布 


Fig.8 Pressure and velocity in inducer outlet 


在 诱导 轮 叶 片 入 口 轮 缘 处 ， 空 化 区 域 狭长 ， 流 体 的 相对 速度 较 大 ， 叶 片 没有 对 流体 做 功 ， 


流体 对 诱导 轮 叶 片 边缘 造成 很 大 的 速度 冲击 ， 使 得 压力 较 低 ， 同 时 诱导 轮 进口 回流 也 会 
造成 局 部 压力 降低 ， 当 压力 低 于 水 的 饱和 蒸汽 压 ， 于 是 会 产生 空 化 区 域 。 

高 速 涡 轮 泵 的 流量 系数 小 ， 容 易 发 生 脱 流失 速 的 问题 。 脱 流失 速 现象 最 先 发 生 在 轮 载 位 
置 处 ， 这 会 导致 在 出 口 处 发 生 压力 和 速度 的 分 布 不 均 的 现象 ， 如 图 8 所 示 ， 又 由 于 诱导 
轮 中 离心 力 的 作用 ， 液 体 被 甩 到 轮 缘 处 ， 所 以 轮 缘 处 压力 大 ， 轮 载 处 压力 小 ， 这 就 会 产 
生出 口 位 置 的 涡流 区 ， 并 且 出 口 涡流 区 和 离心 轮 的 进口 回流 也 有 很 大 关系 。 出 口 涡流 会 


使 得 流体 能 量 损耗 ， 此 部 分 的 总 压 有 所 降低 ， 同 时 出 口 速度 分 布 不 均 ， 在 低压 位 置 的 速 


度 比较 大 ， 于 是 很 容易 发 生 空 化， 使 得 诱导 轮 的 出 口 靠近 轮 载 处 有 一 部 分 空 化 区 域 。 
由 于 诱导 轮 的 出 口 部 分 区 域 有 空 化 ， 并 且 诱 导 轮 的 出 口 存在 涡流 ， 这 都 会 严重 影响 到 离 
心 轮空 化 性 能 ， 因 此 离心 轮 长 叶片 吸力 面 的 进口 处 会 有 空 化 。 


3 空 化 非 定 常 结 


为 了 检测 流 场 压力 脉动 情况 ， 本 文 共 布置 了 20 个 监测 点 ， 如 图 9 所 示 。 监 测 点 1-4 均匀 
分 布 在 诱导 轮 叶 片 进 口 前 , 监测 点 5-8 均匀 分 布 在 诱导 轮 叶 片 出 口 至 离心 轮 入 口 之 间 , 监 
测 点 9-20 均匀 分 布 在 蜗 这 内 部 。 


9 测 点 位 置 定义 


Fig.9 Locations of monitor poi 


监测 点 1-4 均 匀 分 


/ 布 在 诱导 办 进口 


监测 点 5-8 均 匀 分 
布 在 诱导 轮 出 口 


nts 


对 压力 脉动 进行 传 里 叶 变换 (FFT)， 得 到 的 频 域 图 如 图 10 所 示 。 诱 导 轮 前 的 压力 脉动 的 
主 频 为 诱导 轮 叶 片 通 过 频率 。 诱 导 轮 和 离心 轮 之 间 的 压力 脉动 的 主 频 为 诱导 轮 叶片 通过 
频率 ， 并 且 主 频 在 周 向 幅 值 变化 很 大 ， 次 频 为 3 倍 诱导 轮 叶片 通过 频率 ， 同 时 离心 叶片 


通过 频率 也 对 此 区 域 的 月 


E 力 脉动 产生 影响 ， 但 是 变化 不 大 ， 故 在 诱导 轮 叶 片 旋转 对 于 诱 


导 轮 出 口 处 的 压力 脉动 会 造成 很 大 影响 。 蜗 壳 内 各 流 道 出 口 位 置 的 压力 脉动 的 频率 有 诱 
导 轮 、 离 心 轮 叶 片 通过 频率 及 其 二 次 谐 频 ， 从 图 中 发 现 曙 这 内 的 4 倍 转 频 的 幅 值 变化 比 


较 大 ， 测 点 16《〈 位 于 隔 舌 附近 ) 的 4 倍 转 频 的 幅 值 很 高 ， 同 时 离心 轮 叶 片 通过 频率 的 压 


力 脉动 幅 值 变化 比较 大 ， 


这 是 由 于 旋转 的 叶轮 和 蜗 壳 动静 干涉 引起 的 。 
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10 各 测 点 压力 频 域 图 
Fig.10 Pressure at frequency domain 

二 人 公 
4 结论 
采用 ZGB 空 化 模型 可 以 预测 涡轮 泵 内 空 化 流动 。 涡 轮 泵 内 主要 发 生 空 化 的 区 域 为 诱导 轮 
叶片 进口 边缘 处 、 诱 导 轮 出 口 靠近 轮 融 处 、 离 心 轮 叶 片 进口 处 。 本 文 计算 的 涡轮 泵 仍 有 
改进 空间 ， 可 以 通过 改变 诱导 轮 进出 口角 及 离心 叶片 进口 角 来 优化 。 涡 轮 泵 诱导 轮 前 压 
力 脉动 主 频 为 诱导 轮 叶片 通过 频率 ， 诱 导 轮 和 离心 轮 之 间 的 压力 脉动 受 诱导 轮 旋 转 影 响 
较 大 ， 蜗 壳 内 压力 脉动 受 动静 干涉 影响 较 大 。 
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